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(57) Abstract: A DSL modem wilh a receiver (40) and 
a transmitter (30) includes an LPDC (54) encoder which 
utilizes a generation matrix G which is derived from a 
suhstantially deterministic II matrix in order to generate 
redundant parity bits for a block of bits. The II matrix is 
determined by assigning "ones" ol" a first Nj randomly or 
dctcrministically. Then, beginning wilh the second column, 
assignment of "ones" is carried out dctcrministically wilh 
each T in a previous ancestor column generating a "I" 
in the next descendant column based on the rule that 
a descendant is placed one position below or above an 
ancestor. As a result, descending or ascending diagonals 
arc generated. When distributing "ones" in any given 
column, care is taken lo ensure that no rectangles are 
generated in conjunction wilh other "ones" in the current 
column and previous columns. By avoiding generation of 
rectangles, diagonals are interrupted. 
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DSL MODEM UTILIZING LOW DENSITY PARITY CHECK CODES 
BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates generally to telecommunications. More particularly, the 
present invention relates to DSL modems utili2ing low density parity check (LDPQ codes and 
methods of simply generating such LDPC codes. 

2. State of the Art 

LDPC codes were invented by R. Gallager in 1963. R.G. Gallager, "Low-Density- 
Parity-Check Codes", MIT Press, Cambridge, MA 1963. Over thirty years later, a number of 
researchers showed that LDPC code is a constructive code which allows a system to approach 
the Shannon limit. See, e.g., DJ.C. MacKay and R.M. Neal, "Near Shannon limit performance 
of LDPC codes", Electron. Letters, Vol. 32, No. 18, Aug. 1996; DJ.C. MacKay, "Good Error- 
Correcting Codes Based on Very Sparse Matrices", IEEE Transactions on Information Theory, 
Vol. 45, No. 2, March 1999; DJ.C. MacKay, Simon T. Wilson, and Matthew C. Davey, 
"Comparison of Constructions of Irregular Gallager Codes", IEEE Transactions on 
Communications, Vol 47, No. 10, Oct. 1999; Marc P.C. Fossorier, Miodrag Michaljevic, and 
Hideki Imai, "Reduced Complexity Iterative Decoding of LDPC Codes Based on Belief 
Propagation", IEEE Transactions on Communications, Vol. 47, No. 5, May 1999; E. 
Eleftheriou, T. Mittelholzer, and A. Dholakia, Reduced-complexity decoding algorithm for 
LDPC codes", Electron. Letter, Vol. 37, January 2001. Indeed, these researchers have proved 
that LDPC code provides the same performance as Turbo-code and provides a range of trade- 
offs between performance and decoding complexity. As a result, several companies have 
suggested that LDPC code be used as part of the GXite.bis and G.dmtbis standards. IBM 
Corp., "LDPC codes for G.dmtbis and GJithte", ITU - Telecommunication Standardization 
Sector, Document CF-060, Clearwater, Florida, 8-12 Jan. 2001; Aware, Inc., "LDPC Codes for 
ADSL", ITU - Telecommunication Standardization Sector, Document BI-068, Bangalore, India, 
23-27, October 2000; IBM Corp., "LDPC codes for DSL transmission", ITU - 
Telecommunication Standardization Sector, Document BI-095, Bangalore, India, 23-27, October 
2000; IBM Corp., "LDPC coding proposal for G.dmt.bis and G.lite.bis", ITU - 
Telecommunication Standardization Sector, Document CF-061, Clearwater, Florida, 8-12 Jan. 
2001; IBM Corp., Globespan, "G.gen: G.dmtbis: G.Lite.bis: Reduced-complexity decoding 
algorithm for LDPC codes", ITU - Telecommunication Standardization Sector, Document IC- 
071, Irvine, California, 9-13 April, 2001. 

LDPC code is determined by its check matrix H. Matrix H is used in a transmitter 
(encoder) for code words generation and in a receiver (decoder) for decoding the received code 
block. The matrix consists of binary digits 0 and 1 and has size M k *Mj, where M k is the 
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number of columns, and Mj is the number of rows. Each row in the matrix defines one of the 
check equations. If a "1" is located in the k-th column of the j'th row, it means that the k-th bit 
of the code block participates in the j-th check equation. 

Matrix H is a "sparse" matrix in that it does not have many "ones". Generally, the 
matrix contains a fixed number of "ones" Nj in each column and a fixed number of "ones" N k 
in each row. In this case, design parameters should preferably satisfy the equation: M k * Nj 
= Mj * N k (1) 

Although it is convenient to have equal numbers of "ones" in each column and in each row, this 
is not an absolute requirement Some variations of design parameters N k and Nj are 
permissible; i.e., N k (j) and Nj(k) can be functions of j and k, correspondingly. 
In addition, another important constraint for matrix design is that the matrix should not contain 
any rectangles with "ones" in the vertices. This property is sometimes called "elimination of 
cycles with length 4" or "4-cycle elimination". For purposes herein, it will also be called 
Rectangle elimination". 

Generally, there are two approaches in the prior art to designing H matrices. The first 
approach was that proposed by Gallager in his previously cited seminal work, R.G. Gallager, 
"Low-Density-Parity-Check Codes", MTT Press, Cambridge, MA 1963, and consists of a 
random distribution of Nj ones within each matrix column. This random distribution is canied 
out column by column, and each step is accompanied by rectangle elimination within the current 
column relative to the previous columns. The second approach to H-matrix design is based on a 
deterministic procedure. For example, in the previously cited IBM Corp., "LDPC codes for 
G.dmtbis and G.litbis", ITU - Telecommunication Standardization Sector, Document CF-060, 
Clearwater, Florida, 8-12 Jan. 2001, a deterministic H-matrix construction is proposed which 
includes identity matrices and powers of an initial square permutation matrix. 

Both of the prior art approaches to designing H matrices have undesirable characteristics 
with respect to their implementation in DSL standards. In particular, the random distribution 
approach of Gallager is not reproducible (as it is random), and thus, the H matrix used by the 
transmitting modem must be conveyed to the receiving modem. Because the H matrix is 
typically a very large matrix, the transfer of this information is undesirable. On the other hand, 
while the deterministic matrix of IBM is reproducible, it is extremely complex and difficult to 
generate. Thus, considerable processing power must be dedicated to generating such a matrix, 
thereby adding complexity and cost to the DSL modem. Besides, this approach does not allow 
constructing a matrix with arbitrary design parameters M k and Mj. 
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SUMMARY OF THE INVENTION 

It is therefore an object of the invention to provide simple methods of generating 
reproducible H matrices. 

It is another object of the invention to provide DSL modems which utilize simply 
generated reproducible H matrices. 

In accord with these objects which will be discussed in detail below, the DSL modem of 
the invention generally includes a receiver and a transmitter with the transmitter including a 
substantially deterministic LDPC encoder. The encoder is a function of a substantially 
deterministic H matrix (H = A|B) which is determined according to the steps and rules set forth 
below. More particularly, the encoder takes a block of bits and utilizes a generation matrix G = 
A^B which is derived from (i.e., is a function of) the H matrix in order to generate redundant 
parity bits. The redundant bits are appended to the original block of bits to generate a word. 

The substantially deterministic H matrix is determined as follows. First, the "ones" of a 
first column Nj are assigned randomly or deterministically. Preferably, the ones are distributed 
evenly within the first column with the first "1" in the first row of the first column according to 
the algorithm: 

H(r,l)=l, where r=l+(i-l)*integer (Mj/Nj); i=l,2,...Nj (2) 
Then, beginning with the second column, assignment of "ones" is carried out deterministically 
with each "1" in a previous (ancestor) column generating a "1" in the next (descendant) 
column based on the rule that a descendant is placed one position below or one position above 
an ancestor (it being determined in advance by convention whether the position below is used or 
the position above is used). As a result, a descending diagonal or an ascending diagonal is 
generated. Where a descending diagonal is used and the ancestor is in the lowest row of the 
matrix, the descendant may take any position in the next column, although it is preferable to 
place the descendant in the highest free position. 

When distributing "ones" in any given column, each new descendant should be checked 
to ensure that no rectangles are generated in conjunction with other "ones" in the current 
column and previous columns. If a rectangle is generated, the location of the descendant is 
changed, preferably by shifting the location down or up (by convention) one position at a time 
until the descendant is in a position where no rectangle is generated. If the position is shifted 
down and the lowest position is reached without finding a suitable position, the search is 
continued by shifting the location one position up from the initial descendant position until a 
suitable position is found. 
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According to the invention, the descendants may be generated in any given order. Two 
preferable generation orders correspond to increasing or decreasing ancestor positions in the 
column. For example, descendants may be generated by first generating a descendant for the 
ancestor at the bottom of the matrix, then by generating a descendant for the ancestor above that 
in the column, then by generating a descendant for the ancestor above that one, etc. (also called 
herein "bottom-up"); or by first generating a descendent for the ancestor at the top of the 
matrix, then by generating a descendant for the ancestor below that in the column, then by 
generating a descendant for the ancestor below that one, etc. (also called herein "top-down"). 

When generating descendants it is possible that one or more descendants can 
"disappear" because of the lack of free positions satisfying the rectangle elimination criterium. 
To regenerate the "lost descendant", it is generally sufficient to change the order of descendant 
generation for that column. Thus, if the order of descendant generation was conducted "bottom- 
up", the direction of generation is switched to "top-down" and vice versa; preferably for that 
column only. If changing the order of descendant generation in a column does not cause a free 
position to appear, the descendant disappears for that column. 

When a descendant disappears it is desirable in the next column to provide a new 
descendant which does not have an ancestor. In this case, a search of an acceptable position for 
an "ancestor-less" descendant is conducted, preferably from the first row down. 

Additional objects and advantages of the invention will become apparent to those skilled 
in the art upon reference to the detailed description taken in conjunction with the provided 
figures. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a high level block diagram of a DSL modem utilizing LDPC encoding and 
decoding according to the invention. 

Figure 2 is a high level flow diagram of a manner of using an H matrix in the DSL 
modem of Figure 1. 

Figure 3 is a flow chart of a preferred method of generating an H matrix according to the 
invention. 

Figure 4a is an H matrix of size 20x15 generated using bottom-up descendant 
generation. 

Figure 4b is an H matrix of size 20x15 generated using top-down descendant generation. 

Figure 5 is an H matrix of size 276x69 generated using bottom-up descendant 
generation. 

Figure 6 is an H matrix of size 529x69 generated using bottom-up descendant 
generation 

Figures 7a-7c are examples of initialization values for the H-matrix according to a 
second embodiment of the invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Turning to Figure 1, a high level block diagram of a DSL modem 10 is seen. The 
modem 10 preferably includes a digital interface 20, a transmitter section 30 and a receiver 
section 40. The transmitter section preferably includes a scrambler 52 which receives data from 
the digital interface 20, a LDPC encoder 54, an interleave 56, a bit distributor 58, a mapper 60, a 
gain element 62, an infinite fast Fourier transform block (IFFT) 64, a cyclic extension block 66, 
a digital to analog converter 68 and a front end transmit block 69 which interfaces with a hybrid 
70. The receiver section preferably includes a front end receive block 71 which interfaces with 
the hybrid 70, an analog to digital converter 72, a time equalizer (TEQ) 73, a fast Fourier 
transform block (FFT) 74, a frequency equalizer (FEQ) 76, a demapper 78, a deinterleaver 80, a 
LDPC decoder 82, and a descrambler 84 which provides data to the digital interface 20. Other 
than the details of the LDPC encoder 54 (and decoder 82), the modem 10 is substantially as 
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would be understood by those skilled in the art In addition, it will be appreciated by those 
skilled in the art that the modem 10 may be implemented in hardware, software, or a combination 
thereof. 

High level details of the LDPC coder 54 and decoder 82 are seen in Fig. 2. In particular, 
the LDPC coder 54 and decoder 82 utilize an H matrix which is designed according to the steps 
and rules set forth below. The H matrix, where H = A|B, with A being a square matrix and B 
being the remaining matrix rectangle, is used for encoding purposes to generate a generation 
matrix G. Matrix G is defined by G = A^B, which results from multiplying the inverse of the 
square A matrix with the rectangular B matrix. The LDPC encoder 54 uses the G matrix and a 
block of bits received from the scrambler 52 to generate a set of parity bits (also called redundant 
bits). The parity bits are appended to the block of bits received from the scrambler 52 to 
generate a word which is forwarded to the interleaver 54 and further processed. If desired, and 
as suggested by Fig. 2, rather than appending the redundant bits to the block of data, the G 
matrix may include an "identity matrix" portion so that a multiplication of the G matrix and the 
block of bits directly provides the resulting word 

The H matrix is likewise used on the decoding side. In particular, deinterleaved words 
received by the LDPC decoder are subjected to soft decisions (as is known in the art), and then 
subjected to probabilistic decoding which requires information of the H matrix which was 
utilized to generate the parity bits. 

The H matrix (and G matrix) may be generated by a microprocessor (not shown) and 
software which may also be used to implement one or more additional elements of the 
transmitter or receiver of the modem 10. Alternatively, the H matrix (and G matrix) may be 
implemented in other hardware and/or software in the modem 10. Technically, only the G matrix 
needs to be available for the transmitter (encoding), while only the H matrix is needed for the 
receiver (decoding). 

According to the invention, the H matrix is a substantially deterministic matrix which, 
according to a first embodiment, may be determined according to the steps of Fig 3. First, at 
step 102 the "ones" of a first column Nj are assigned randomly or deterministically. 

Preferably, the "ones" are distributed evenly within the first column with the first "1" in the 
first row of the first column according to relationship (2) set forth above: 

H(r,l)=l, where r=l+(i-l)*integer (Mj/Nj); i=l,2,...Nj 
where Mj is the number of rows in the matrix and Nj is the number of "ones" in the column. 
Thus, if the "ones" are assigned deterministically, the first "one" is located at H(l,l) and the 
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remainder of "ones'* for the column are evenly distributed in the column. If, on die other hand, 
the "ones" are assigned randomly, preferably, a "one" is located in a random row of column 1, 
and the remaining "ones" are evenly distributed While less preferred, all "ones" in column 
one can be randomly located. 

Returning to Figure 3, once the "ones" of the first column are assigned, at 103 the next 
column is addressed. In particular, at 104, each of the "ones" of the next column is generated 
deterministically (i.e., according to a predetermined set of rules). In particular, a "one" of the 
second column (called a "descendant") is generated at 104 by placing the descendant "1" one 
position below or one position above its "ancestor" "one" of the previous column (it being 
determined in advance by convention whether the position below is used or the position above is 
used). As a result, a descending diagonal or an ascending diagonal is generated. Where a 
descending diagonal is used and the ancestor is in the lowest row of the matrix, the descendant 
may take any position in the next column, although it is preferable to place the descendant in the 
highest free position. This may be seen with reference to columns 5 and 6 of Figure 4b. As 
seen in Fig. 4b, H(15,5) = 1, and accordingly, the descendant is found in the first row of column 
6; i.e., H(l,6) = 1. Similarly H(15,9) generates H(l,10), and H(15,12) generates H(l,13). 
Conversely, where an ascending diagonal is used and the ancestor is in the highest row of the 
matrix, the descendant may take may position in the next column, although it is preferable to 
place the descendant in the lowest free position. 

When distributing "ones" in any given column, at 106, each new descendant is checked 
to ensure that no rectangles are generated in conjunction with other "ones" in the current 
column and previous columns. If a rectangle is generated, a command to change the location of 
the descendant is issued at 108, preferably by shifting the location down or up (by convention) 
one position at a time (at 104) until the descendant is in a position where no rectangle is 
generated (as determined at 106). If the position is shifted down and the lowest position is 
reached without finding a suitable position, the search is continued by shifting the location one 
position up from the initial descendant position until a suitable position is found. 

Rectangle elimination is seen in the matrix of Figure 4a, In particular, referring to the 
fifth and sixth columns, according to the rule of descendants, ancestor H(5,5) = 1 should 
generate a descendant H(6,6) = 1. However, this descendant would cause a rectangle to appear 
in conjunction with H(l,6), H(l,l), and H(6,l). Going down in column 6, it is seen that position 
H(7,6) is also not acceptable as it would cause a rectangle to appear in conjunction with H(7,2), 
H(12,2) and H(12,6). Thus, the descendant of H(5,5) is found in position H(8,6). 
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According to the invention, the descendants may be generated in any given order. Two 
preferable generation orders correspond to increasing or decreasing ancestor positions in the 
column. For example, descendants may be generated by first generating a descendant for the 
ancestor at the bottom of the matrix, then by generating a descendant for the ancestor above that 
in the column, then by generating a descendant for the ancestor above that one, etc. (also called 
herein "bottom-up"). The bottom-up technique is seen in Fig. 4a, where a descendant is 
generated first for H(15,5), then for ancestor H(10,5), and finally for ancestor H(5,5). 
Alternatively, descendants may generated by first generating a descendent for the ancestor at the 
top of the matrix, then by generating a descendant for the ancestor below that in the column, then 
by generating a descendant for the ancestor below that one, etc. (also called herein "top-down"). 
The top-down technique generates a fiill diagonal of "ones" from H(l,l) to H(Mj,Mj) as is 

seen in Fig. 4b. In Fig. 4b, a descendant is generated first for H(5,5), then for ancestor H(10,5), 
and finally for ancestor H(15,5). Regardless of whether the top-down or bottom-up technique is 
used, at 110, a determination is made as to whether all descendants for a column have assigned. 
If all descendants for the column have not been assigned, the program cycles through steps 104- 
110. If all descendants for the column have been assigned, unless a determination is made at 
1 12 that the column is the last column, the next column is accessed at 103 for placement of 
descendant "ones". 

When generating descendants it is possible that one or more descendants can 
"disappear*' because of the lack of free positions satisfying the rectangle elimination criterium. 
This determination can be made at step 1 15 (shown in phantom after step 108). To regenerate 
die "lost descendant", it is generally sufficient to change the order of descendant generation for 
that column (step 1 17 - shown in phantom). Thus, if the order of descendant generation was 
conducted "bottom-up", the direction of generation is switched to "top-down" and vice versa. 
Preferably, the order of descendant generation is changed only for that column. If changing the 
order of descendant generation in a column does not cause a free position to appear, the 
descendant disappears for that column. 

When one or more descendants disappear in a column, it is desirable in the next column 
to provide a new descendant for each descendant which does not have an ancestor. In this case, a 
search of acceptable positions for each "ancestor-less" descendant is conducted, preferably 
from the first row down. 

Generally, as set forth above, the number of "ones" in each column Nj is determined by 
the number of "ones" in the previous column, because a descendant is generated for each 
ancestor. In the preferred embodiment of the invention, this number is fixed and defined as a 
design parameter. On the other hand, the number of "ones" in each row (row weight) is 
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preferably limited to a maximum row weight (Max(N k )) which is also a design parameter. Thus, 
if during distribution of "ones" within a particular column, the number of "ones" in some row 
reaches Max(N k ), "ones" should not be inserted in that row (i.e., the remaining part of the row 
is automatically filled with "zeros"), and the descendant "one" is moved by shifting the 
location of the descendant one position down or up (by convention). 

An implementation in Matlab of the method of H matrix design according to Fig. 3 as 
described above is as follows: 
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Check-Matrix initialization 



Mkoinput ( 1 number of matrix columns, code block length 

NjoinputP number of "ones* 1 in a column, number of checks for bit 

Nk» input ('number of "ones" in a row, number of bits in each check 

Mj=input( •number of matrix rows, number of check equations 



Mk» 


'); 


Nj* 


•); 


Hk= 


•>; 


Mj« 


') ; 



c«l); 

w«[0 2*(ones(9ize(lt(Mk-l)>))]; 
for j«l:Mj 
C«[C;w] ; 

end 

vNkozeros (size (l:Mj) ) ; 
vNjozeros (size(l:Mk) ) ; 

1-st column initialization 
rr=floor (Mj/Nj) ; 
for jr«l:Nj 

r«l+ (jr-l) *rr; 

C(r,l)«l; 

vNk(r)^l; 

end 

vNj (l)«Nj; 

Matrix Design 
for k=l: (Mk-1) 
Z=C(: , (k+1) ) ; 
for h*l:2 
if h-»l 

count=0 / counthl«0 ; 
for jj«=l:Mj 
X»0 ; 

j*Mj+l-jj; 

if j— Mj Sc C(Mj,k)==l 
n=3 0;nn»0; 

while nrinoO & n<Mj 
n=n+l ; 

if C(n, (k+l) )=«2 

C(n, (k+l))«l; 

nn«l ; counthl=counthl+l ; 

end 

end 

xen; 

elseif C(j,k)«=l & C((j+1) , {k+l))*«2 

C((j+1) . (k+1) 

x=j+l; 

counthl«counthl+l ; 

elseif C(j,k)«i & c( <j+l) , (k+1) ) <2 
m«0;nn=0; 

while nn««0 & ra<(Mj-l) 
m»nw-l ; 

ifi (j+l+m) < (Mj+1) 
nm»m; 

elseif (j+X+m)>Mj 
nm=Mj -j -l-m; 



%Check-Matrix 



%current numbers of "ones" in rows 
%current numbers of "ones** in columns 



% evenly distributed ones 



%column by column "1" assignment 

%h=*l: search, beginning from the last row 

Icurrenc number of "ones" in the column 
%row by row assignment, beginning from Mj 

% transfer "l" from last row to 1-st row 



%typical diagonal shift 

%additional shift 

%searching the acceptable place 

% searching down 

% searching up 
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end 

if C ( j+i+nm, (k*i) ) «o2 

C(j+i+nm, (k+i))-i; 

nn= 1 ; count hi ocounthl+1 ; 

end 

end 

xmj+l+nm; 
end 

if x>0 %rectangle elimination 

count =*count+l; 

kkok; 

while kk>0 
if C(x,kk) «nl 
for jjolrMj 

if <C(jj,kk)«X) & (abs<jj-x)>0) & (count<Nj) 

C(jj, (k+l))n:0; 

end 

end 

end 

kk«=kk-l; 

end 

end 

end %for j j~l :Mj (end of one column design for h»l) 

elseif h«2 & counthl<Nj 
count «0 ; counth2»0 ; 

for jj»l:Mj %row by row "l n assignment from l-st row 

x=0; 

j-jj; 

if j«»Mj & C(Mj,k)««l %transfer *1" from last row to l-st row 

n«0 ; nn=0 ; 

while nna«0 & n<Mj 

nan+1 ; 

if z(n)==2 

Z(n)**l; 

nn«l ; 

count h2 scoun th2 +1 ; 

end 

end 

x=n; 

elseif C(j,k)»»l & a(j+l)«=2 

counth2=counth2+l ; 

elseif C(j,k)«=l & z{j+l)<2 

m«0 ;nn=0 ; 

while nn»=0 6 m<(Mj-l) % searching the acceptable place 

oiBfii+l; 

if (j+l+m) <(Mj+l) % searching down 

IUR°R1; 

elseif (j+l+m)>Mj %searching up 

nm«Mj - j -1-ra; 

end 

if z(j+l+nra)-«-2 
2(j+l+nm) -1; 
nn»i; 

counth2»counth2+l j 
end 
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%rectangle elimination 



end 

x*j+l+nm; 
end 
if x>0 

count = count +1 t 
kkok; 

while kk>0 
if C(x,kk)o=l 
for jj«l:Mj 

if (C{jj,kk)»ol) & (abs(jj-x)>0) & (count<Nj) 

z(jj)=0; 

end 

end 

end 

kk=kk-l; 
end 
end 
end 

if counth2>counthl 
C(: , (k+1) )«z; 
end 
end 
end 

if vNj (k)<Nj 
qq«0;f =0; 
while f<l & qq<Mj 

if C(qq. (k+l))~2 
C<qq, (k+l))«l; 
f~f+l 
end 
end 
end 

for jj*l:Nj 
if C(jj. (k+l)>«l 
vNk(jj)«vNk(jj)+l; 
vNj (k+ 1) =vNj (k+1) +1; 
else 

C(jj, (k+m-0; 
end 

if vNk(jj)»»Nk 
for xk»(k+2) :Mk 
C(jj,kk>»0? 
end 
end 
end 
end 
C; 
vNj 
vNk 



%for jj«i:Mj(end of one column design for h=2 



%if h==l 
%for h=l:2 

%ancestor recreation 



%calculation of ones in each row 
%calculation of ones in each column 

%change w 2" to "0 W in columns 



%change W 2 W to "0 1 in rows 



%for k°l( (Mk-1) (end of columns design) 
Idemot Check Matrix 
tdemo: Number of ones in columns 
% demo: Number of ones in rows 
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It will be appreciated by those skilled in the art that other implementations of generating 
an H matrix design in Madab or in other software or hardware are easily obtained. 

Figure 5 is an H matrix of size 276x69 generated using bottom-up descendant 
generation as set forth in the previously listed Madab program. The H matrix of Figure 5 has 
design parameters M k = 276, Mj = 69, N k = 12, Nj = 3. The generated H matrix contains a 
fixed number of "ones" in columns Nj = 3, and a fixed number of "ones" in rows N k = 12. 
Similarly, Figure 6 is an H matrix of size 529x69 generated using the previously listed Madab 
program. The H matrix of Figure 6 has the design parameters M k = 529, Mj = 69, Max(N k ) = 
12, Nj = 3. This H matrix contains a fixed number of "ones" in its columns, Nj = 3; and a 
nonfixed, but limited number of "ones" in its rows, 20 < N k < 25. In Figures 5 and 6, a "one" 
is shown by a dot, while a "zero" is shown by the absence of a dot. 

According to a second embodiment of the invention, the design procedure for generating 
the H matrix may be simplified. In particular, because every column should contain at least one 
"1", it is possible to initialize the H matrix with an effectively continuous diagonal. Three such 
diagonals are shown in Figures 7a-7c, with Fig. 7a representing a descending diagonal, Fig. 7b 
representing an ascending diagonal, and Fig. 7c representing a mixed descending-ascending 
diagonal. Of course, an ascending-descending diagonal (not shown) could likewise be utilized. 
With the H matrix initialized as shown, the steps shown in Fig. 3 are carried out only with 
respect to the "ones" which are distributed in the first column and their descendants, thereby 
reducing the number of calculations required. 

With the substantially deteraunistic method of generating H matrices set forth above, it 
will be appreciated that if standard conventions (e.g., detenninistic first column, descending 
diagonal generation, bottom-up descendant generation) are agreed upon for all modems, the only 
information which must be transferred from a transmitting modem to a receiving modem 
regarding the H matrix includes the matrix size (M k x Mj), and the number (or maximum 
thereof) of "ones" in a row or column; N k and Nj. If standard conventions are not used, code 
representing one or more of: whether descending or ascending diagonals are used, whether 
bottom-up or top-down descendant generation is used, the basis of the first column, etc. will also 
be required to be sent from the transmitting modem to the receiving modem Regardless, the 
generation of the H matrix (and hence the G matrix) will be gready simplified in both the 
transmitter and receiver. 

There have been described and illustrated herein embodiments of modems utilizing 
LDPC coders based on particular H matrices, and methods of simply generating such H 
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matrices. While particular embodiments of the invention have been described, it is not intended 
that the invention be limited thereto, as it is intended that the invention be as broad in scope as the 
art will allow and that the specification be read likewise. Thus, while particular code has been 
listed for generating H matrices, it will be appreciated that other software and/or hardware could 
be utilized. Also, while the H matrix was discussed with reference to a particular DSL- type 
modem, it will be appreciated that the H matrix could be used in other types of modem or in 
other applications. Further, while particular preferred conventions were described, it will be 
appreciated that other conventions could be added or substituted. For example, while a "bottom- 
up** and a "top-down" convention were described, a "middle-out" convention could be 
utilized. Similarly, while the convention of causing the descendant to be located in a row one 
position down or up from the ancestor of the previous column is preferred, a diagonal can be 
likewise generated by causing the descendant to be located two, three or n rows up or down from 
the ancestor of the previous column. In addition, the convention utilized to generate the 
descendants could change, by convention, from column to column. Furthermore, while rectangle 
elimination is shown in Fig. 3 to be conducted upon placement of each "1" value in the matrix, 
it will be appreciated that no checking is required for the first few columns which in principle 
cannot create a rectangle. Also, while Fig. 3 represents checking for rectangle elimination after 
each placement of a "1", it is equivalently possible (and is in fact shown in the Matlab program 
described above) to determine in advance for each column, into which rows a "1" value cannot 
be placed due to the rectangle rule. Thus, many equivalent flow charts such as Fig. 3 may be 
generated which represent methods of generating an H matrix according to the invention. 
Further yet, while the invention was described as generating a matrix by inserting "1" values 
into a first column of the matrix and assigning descendant ones in subsequent columns, it will be 
appreciated that the "1" values could be inserted from left to right, or from right to left in the 
matrix, and the first column to received the ones could be any column of the matrix. Where, a 
middle column is selected as the first column to receive the ones, the first and last columns will 
be perceived to be adjacent each other for purposes of continuing the assignment of descendant 
ones. It will therefore be appreciated by those skilled in the art that yet other modifications could 
be made to the provided invention without deviating from its spirit and scope as so claimed. 
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1. A digital subscriber line (DSL) modem, comprising: 

a) a digital interface; and 

b) a transmitter coupled to said digital interface, said transmitter including a low density parity 
check (LDPC) encoder which generates redundant bits utilizing a generation matrix which is a 
function of a substantially deterministically generated H matrix; 

c) a receiver coupled to said digital interface, said receiver including a LDPC decoder; and 

d) means for substantially deterministically generating said H matrix, said H matrix having a 
plurality of columns (M k ) and a plurality of rows (Mj), said means for generating said H matrix 
being associated with at least one of said transmitter and said receiver and including means for 
assigning a plurality of "ones" into a first column of said H matrix, means for deterministically 
assigning descendant "ones" in adjacent columns in a substantially diagonal fashion based on 
the assignment of ancestor "ones" in said first column, and means for preventing "ones" from 
being located in locations which would cause a rectangle pattern of "ones" in said H matrix by 
causing a descendant "one" which would otherwise be located in one of said locations to 
assume another location, thereby causing breaks in diagonal patterns of "ones". 

2. A modem according to claim 1, wherein: 

said first column of the H matrix is a left-most column of said matrix. 

3. A modem according to claim 1, wherein: 

said adjacent columns are subsequent columns. 

4. A modem according to claim 1, wherein: 

said means for generating an H matrix assigns "ones" into said first column in a 
random fashion. 

5. A modem according to claim 1, wherein: 

said means for generating an H matrix assigns "ones" into said first column in a 
deterministic fashion, 

6. A modem according to claim 1, wherein: 

said means for generating an H matrix assigns "ones" into said first column such that 
said "ones" are distributed evenly within said first column. 
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7. A modem according to claim 1, wherein: 

said means for generating an H matrix assigns "ones" into said first column according 

to 

H(r,l)=l, where r=l+(i-l)*integer (M/Nj); i=l,2,...Nj 

where Mj is the number of rows in said H matrix, and Nj is the number of "ones" in a column 
of said H matrix. 

8. A modem according to claim 1, wherein: 

said means for deterministically assigning descendant "ones" in adjacent columns 
locates a descendant one position below an ancestor unless locating said descendant in that 
position which would cause a rectangle pattern of "ones" in said H matrix, 

9. A modem according to claim 8, wherein: 

when an ancestor in a particular column is located at a last row of said H matrix, said 
means for deterministically assigning descendant "ones" locates a descendant for that ancestor 
in any row of the next column which would not cause a rectangle pattern of "ones" in said H 
matrix. 

10. A modem according to claim 9, wherein: 

said any row of the next column is the highest row possible without causing a rectangle 
pattern of "ones" in said H matrix. 

11. A modem according to claim 1, wherein: 

said means for deterministically assigning descendant "ones" in adjacent columns 
locates a descendant one position above an ancestor unless locating said descendant in that 
position which would cause a rectangle pattern of "ones" in said H matrix. 

12. A modem according to claim 1 1, wherein: 

when an ancestor in a particular column is located at a first row of said H matrix, said 
means for deterministically assigning descendant "ones" locates a descendant for that ancestor 
in any row of the next column which would not cause a rectangle pattern of "ones" in said H 
matrix. 

13. A modem according to claim 1, wherein: 

said means for preventing "ones" from being located in locations which would cause a 
rectangle pattern of "ones" in said H matrix comprises means for shifting the location of a 
"one" by one position at a time until the descendant is in a position where no rectangle is 
generated. 
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14. A modem according to claim 13, wherein: 

said means for preventing shifts the location down by one position at a time. 

15. A modem according to claim 14, wherein: 

if a lowest row is reached and still does not present a suitable location for the descendant, 
said means for preventing shifts the location up by one position at a time. 

16. A modem according to claim 13, wherein: 

said means for preventing shifts the location up by one position at a time, 

17. A modem according to claim 16, wherein: 

if a highest row is reached and still does not present a suitable location for the 
descendant, said means for preventing shifts the location down by one position at a time. 

18. A modem according to claim 1, wherein: 

said means for deterministically assigning descendant "ones" in adjacent columns in a 
substantially diagonal fashion based on the assignment of ancestor "ones" in said first column 
generates descendants in a bottom-up manner such that said means for deterministically 
assigning descendant "ones" first generates a descendant for a particular ancestor located in the 
lowest row of said first column containing an ancestor, and then generates a descendant for the 
ancestor above that particular ancestor. 

19. A modem according to claim 18, wherein: 

said means for deterministically assigning descendant "ones" switches from generating 
descendants in said bottom-up manner to generating descendants in a top-down manner if said 
means for preventing would otherwise prevent a descendant from being located in any location 
due to a rectangle pattern of "ones" appearing in said H matrix. 

20. A modem according to claim 19, wherein: 

said means for deterministically assigning descendant "ones" 
generates an ancestor-less descendant in a column adjacent a column where a descendant 
disappears because it cannot be located by said means for deterministically assigning descendant 
"ones" because of said means for preventing. 
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21. A modem according to claim 1, wherein: 

said means for deterministically assigning descendant "ones" in adjacent columns in a 
substantially diagonal fashion based on the assignment of ancestor "ones" in said first column 
generates descendants in a top-down manner such that said means for deterministically 
assigning descendant "ones** first generates a descendant for a particular ancestor located in the 
highest row of said first column containing an ancestor, and then generates a descendant for the 
ancestor below that particular ancestor. 

22. A modem according to claim 21, wherein: 

said means for deterministically assigning descendant "ones" switches from generating 
descendants in said top-down manner to generating descendants in a bottom-up manner if said 
means for preventing would otherwise prevent a descendant from being located in any location 
due to a rectangle pattern of "ones" appearing in said H matrix, 

23. A modem according to claim 22, wherein: 

said means for deterministically assigning descendant "ones" 
generates an ancestor-less descendant in a column adjacent a column where a descendant 
disappears because it cannot be located by said means for deterministically assigning descendant 
"ones" because of said means for preventing, 

24. A modem according to claim 1, wherein: 

said transmitter further includes a scrambler coupled to said digital interface and to said 
LDPC encoder 54, a mapper coupled to said LDPC encoder, a digital to analog converter 
coupled to said mapper, and a hybrid coupled to said digital to analog converter. 

25. A modem according to claim 24, further comprising: 

said receiver includes a demapper coupled to said hybrid, wherein said LDPC decoder is 
coupled to said demapper, 

26. A modem according to claim 24, wherein: 

said H matrix is defined by H = AjB, with A being a square matrix and B being a 
remaining matrix rectangle, and 

said generation matrix is defined by G = A^B, where A 4 is an inverse of said square 
matrix A. 

27. A modem according to claim 26, wherein: 

said LDPC encoder uses said generation matrix and a block of bits received from said 
scrambler to generate a set of parity bits. 



BNSOOCID: <WO 02095985A1 . 1_> 



WO 02/095965 



PCT/US02/11835 



19 

28. A method of generating an H matrix for a low density parity check code, comprising: 

a) assigning a plurality of "ones" into a first column of an H matrix having a plurality of 
columns (M k ) and a plurality of rows (Mj); 

b) based on the assignment of "ones" in said first column, deterministically assigning 
descendant "ones" in adjacent columns in a substantially diagonal fashion; 

c) preventing "ones" from being located in locations which would cause a rectangle pattern of 
"ones" in said H matrix by causing a descendant "one" which would otherwise be located in 
one of said locations to assume another location, said preventing thereby causing breaks in 
diagonal patterns of "ones". 

29. A method according to claim 28, wherein: 

said first column of the H matrix is a left-most column of said matrix. 

30. A method according to claim 28, wherein: 

said adjacent columns are subsequent columns. 

31. A method according to claim 28, wherein: 

said assigning "ones" into said first column comprises assigning said "ones" into said 
first column in a random fashion. 

32. A method according to claim 28, wherein: 

said assigning "ones" into said first column comprises assigning said "ones" into said 
first column in a deterministic fashion. 

33. A method according to claim 28, wherein: 

said assigning "ones" into said first column comprises assigning said "ones" such that 
said "ones" are distributed evenly within said first column. 

34. A method according to claim 28, wherein: 

said assigning "ones" into said first column comprises assigning said "ones" 
according to 

H(r,l)=l, where r^l+(i-l)*integer (Mj/Nj); i=l,2,...Nj 
where Mj is the number of rows in said H matrix, and Nj is the number of "ones" in a column 
of said H matrix. 



BNSDOCID: <WO 02095965A1..I_> 



WO 02/095965 



PCT/US02/11835 



20 

35. A method according to claim 28, wherein: 

said assigning descendant "ones" in adjacent columns comprises locating a descendant 
one position below an ancestor unless locating said descendant in that position which would 
cause a rectangle pattern of "ones" in said H matrix. 

36. A method according to claim 35, wherein: 

when an ancestor in a particular column is located at a last row of said H matrix, said 
assigning descendant "ones" comprises locating a descendant for that ancestor in any row of 
the next column which would not cause a rectangle pattern of "ones" in said H matrix. 

37. A method according to claim 36, wherein: 

said any row of the next column is the highest row possible without causing a rectangle 
pattern of "ones" in said H matrix. 

38. A method according to claim 28, wherein: 

said assigning descendant "ones" in adjacent columns comprises locating a descendant 
one position above an ancestor unless locating said descendant in that position which would 
cause a rectangle pattern of "ones" in said H matrix. 

39. A method according to claim 38, wherein: 

when an ancestor in a particular column is located at a first row of said H matrix, said 
assigning descendant "ones" comprises locating a descendant for that ancestor in any row of 
the next column which would not cause a rectangle pattern of "ones" in said H matrix. 

40. A method according to claim 28, wherein: 

said preventing "ones" from being located in locations which would cause a rectangle 
pattern of "ones" in said H matrix comprises shifting the location of a "one" by one position 
at a time until the descendant is in a position where no rectangle is generated. 

41. A method according to claim 40, wherein: 

said preventing comprises shifting the location down by one position at a time. 

42. A method according to claim 41, wherein: 

if a lowest row is reached and still does not present a suitable location for the descendant, 
said preventing comprises shifting the location up by one position at a time. 

43. A method according to claim 40, wherein: 

said preventing comprises shifting the location up by one position at a time. 
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44. A method according to claim 43, wherein: 

if a highest row is reached and still does not present a suitable location for the 
descendant, said preventing comprises shifting the location down by one position at a time. 
]45. A method according to claim 28, wherein: 

said deterministically assigning descendant "ones" in adjacent columns in a 
substantially diagonal fashion comprises generating descendants in a bottom-up manner such 
that a descendant is first generated for a descendant for a particular ancestor located in the lowest 
row of said first column containing an ancestor, and then a descendant is generated for the 
ancestor above that particular ancestor. 

46. A method according to claim 45, wherein: 

said deterministically assigning descendant "ones" switches from generating 
descendants in said bottom-up manner to generating descendants in a top-down manner if said 
preventing would otherwise prevent a descendant from being located in any location due to a 
rectangle pattern of "ones" appearing in said H matrix. 

47- A method according to claim 46, wherein: 

said deterministically assigning descendant "ones" 
generates an ancestor-less descendant in a column adjacent a column where a descendant 
disappears because it cannot be located at all because of said preventing. 

48. A method according to claim 28, wherein: 

said deterministically assigning descendant "ones" in adjacent columns in a 
substantially diagonal fashion comprises generating descendants in a top-down manner such that 
a descendant "one" is first generated for a descendant for a particular ancestor located in the 
highest row of said first column containing an ancestor, and then a descendant is generated for 
the ancestor below that particular ancestor. 

49. A method according to claim 48, wherein: 

said deterministically assigning descendant "ones" switches from generating 
descendants in said top-down manner to generating descendants in a bottom-up manner if said 
preventing would otherwise prevent a descendant from being located in any location due to a 
rectangle pattern of "ones" appearing in said H matrix, 

50. A method according to claim 49, wherein: 

said deterministically assigning descendant "ones" 
generates an ancestor-less descendant in a column adjacent a column where a descendant 
disappears because it cannot be located at all because of said preventing. 
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